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ABSTRACT 



^ I VFTS102 has a projected rotational velocity (>500km s ^) and would appear 



to be the fastest rotating O-type star. We show that its high rotational velocity 



■ could be understood within the framework of the binary merger. In the binary 

merger channel, the progenitor binary of VFTS102 would evolve into contact 



while two components are still on the main sequence, and then merge into 
a rapidly rotating single star. Employing Eggleton's stellar evolution code, 
we performed binary stellar evolution calculations and mapped out the initial 
parameters of the progenitor of VFTS102 in the orbital period-mass ratio (P — 
q) plane. We found that the progenitor binary of VFTS102 with initial mass 
ratio go ^ 0.7 should have an initial orbital period shorter than 3.76—4.25 days, 
while above this mass ratio it should have an initial orbital period shorter than 
1.44 — 1.55 days. The progenitor of VFTS102 would evolve into contact during 
the rapid mass transfer phase or during the subsequent slow mass transfer 
phase, and might ultimately merge into a rapidly rotating massive star. In 
addition, we performed Monte Carlo simulations to investigate the binary 
merger channel. We estimated the fraction of binaries that would merge into 
single stars and the fraction of single stars that might be produced from the 
binary merger channel. It is found that about 8.7% of binaries would evolve 
into contact and merge into rapidly rotating single stars, and about 17.1% of 
single stars might be produced from the binary merger channel and should 
have similar properties to VFTS102. This suggests that the binary merger 
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channel might be one of the main channels for the formation of rapidly rotating 
massive stars like VFTS102. 

Key words: instabilities - stars: early-type - stars: formation - stars: evo- 
lution - stars: rotation 



1 INTRODUCTION 



VFTS102 is a rapidly rotating 0-type star in the 30 Doradus region of the Large Magellanic 
Cloud, and its projected rotationa l velocity is larger than 500 km s~^ and probably as large 



as 600 km s 



Dufton et al 



2011). 



large surveys (IMartayan et al. 



2006 



t rotates more rapidly tha n any observed stars in recent 



most rapidly rotating massive star (IDufton et al. 



Hunter et a 



2008 



20091 ) and would appear to be the 



201ll ). The rapidly rotating stars such as 



VFTS102 form an important class of objects in several respects. These stars could help 
us to st udy the rotational effects that are important, from star formation through to their 



deaths ( Hunter et al.1 



20081 ). Fast rotation could change the lifetime a nd evolution of rnassive 



stars by rotationa. 



Meynet fc Maedei 



stellar surface ( iHunter et al. 



y ind uced mixing between the core and envelope ( iHeger &: Langer 



2000 : 



2000f). and produce enrichment of a number of different elements at th e 



2009 



Frischknecht et al 



2010 



Potter. Tout fc Eldridgell2012[ ). 



'he rotational 



( jCeorgy et al. 



y ind uced mixing were used to explain the variety of core collapse supernovae 



20091 ). In addition, the rotating ma ssive stars might be the progenitors of 



gamma-ray 



2006 



Tout 



bursts through homogeneous evolution (lYoon fc Langer 



2011 



2005 



Woosley fc Heger 



Two channels for the formation o f VFTS102 have be en proposed in the past. The first 



channel is the mass transfer channel fiDufton et al. 



201 ll ). In this channel, VFTS102 might 



be the mass gainer and be spun up by a past episode of Roche lobe overflow (RLOF ) 



in an interac t ing b inary system (jPacket 



1981 



Cantiello et al 



2007; iDufton et al 



Dufton et al. 



201lh . 



( I2OIII ) suggested that VFTS102 became a runaway star after its compan- 



ion exploded as a supernova. This channel could also explain that the radial velocity of 



VFTS102 differs by 40 km s'^ 
channel is the collision channel 



Fujii. Saitoh &: Portegies Zwart 



Tom the mean radial 



iFrver fc Heger 



2005 



velocity for 30 Doradus. The secon d 



Fujii. Saitoh fc Portegies Zwartll2012l ). 



(I2OI2I ) proposed that VFTS102 had experienced a earlier 



collision, and then might be ejected from the cluster center. The fast rotation of VFTS102 
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might be the result of the colhsion with another star in the cluster center. In addition , 



VFTS102 might be born as rapid rotator (IHuang. Gies fc McSwain 



201Q 



BrottetaL 



20111 ) 



or its surface might spin up during the main sequence (MS) ev olution because angu lar mo- 
mentum was transported from the center to the stellar surface (lEkstrom et al.ll2008l ). 

From a binary evolution point of view, these formation channels are not complete. The bi- 
nary merger might be another possible channel to produce the fastest rotating star VFTS102. 
In this channel, the initial detached binary would experience case A mass transfer and evolve 
into contact while two components are still on the MS, and then merge into a rapidly ro 



fatin g single star (jWebbink 



1976 



Jiang et al. 



2010 



de Mink. Langer fc Izzard 



2011 



Langer 



20121 ). It is w ell known that the 



is very large (IMason et al. 



2009; 



raction. of ni a ssive stars that ar e members of a close binary 
Langerll2012l ). ISana et all ^20121 ) found that 20% of all stars 



born as 0-type stars would merge as a result of case A mass transfer. In addition, the ob- 



('P=1.495d 



Qian et al 



Lorenz. Maver. fc Drechsell Il999h and TU Mus (P=1.387d 



, such as V60 


6 Cen 


Terrell et al. 


2003; 



20071), and ^ 5% of 0-type stars appear to be in contact with orbital pe riods of 



1 — 5.5 days ( iGarmany. Conti fc Massey 



1980; 



Hilditch fc Bell 



1987 



Eggletonlll996h . These 



massive contact binaries might merge into single stars in a similar way as low-mass contact 
binaries, which might merge into single stars due to tidal instability when the sp in angular 



momenturn of the system is m ore than a third of its orbital angular momentum ( iHut 



1980; 



Rasio 



1995 



Li fc Zhang! 120061). or due to th ermal instability when the primary attempts to 



cross the Hertzsprung gap (IWebbink 



1976). 



to result in single, massive stars ( IChen &: Han 



he rnerger of co ntact binaries are expected 



2008 



Tout 



201ll ). which wou 



rapid 



2007 



y rotating due to the orbital angular momentum of the binary sy stems ( ICantiello et al. 



Jiang et al. 



2010 



d be extremely 



de Mink. Langer fc Izzard 



2011 



Langer 



20121 ). Therefore, the binary 



merger might be a possible channel for the production of VFTS102. 

The evolution of massive clo s e binaries h as been investigated by many au t hors, e.g. 



Podsiadlowski. Joss fc Hsul (11992 ) :IPols ( 



de Mink. Pols fc Hilditchl (120071 1. 



994) 



Wellstein. Lariger fc BraunI (l200l[ ): lNelson &: Eggleton 



Nelson fc EggletonI (120011 ) calculated the case A bi- 



nary evolution and constructed a large grid of models (0.8< Mio < 50). They mapped out 
the initial parameters in the mass ratio-orbital period plane of six subtypes of case A evolu- 
tion. They suggested that three of these subtypes (AR, AS, AD) lead to contact while both 
components are on the MS. In case AR or AS, the secondary expands in response to the 
thermal time-scale mass transfer or the nuclear time-scale mass transfer from the primary 
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and fills its own Roche lobe. The system probably form a stable contact binary, and might 
ultimately merge into a single star. However, in case AD, the secondary cannot accrete all 
the proffered material transfer from the primary on t he dynamical time-scale. This probably 



leads very quickly to a common envelope, spiral- in (jPaczyhski 



1987 



Ge et al 



Jiang et al. 



2001 



by 



1976 



Hiellming fc Webbink 



20101), and coalescence on a quite short timescale (INelson &: Eggleton 



2012 



2001 



Jiang et al. 



. Case AD might be common in the low mass binaries (INelson fc Eggleton 



Becker et al.l (1201 ll ) and 



2012), espe cially in th o se bin aries with two M dwarf components observed 

(I2OI2I ). Few massive syst ems are classified as AD 



Nefs et al 



and most of massive contact sy stems are classified as AR or AS (INelson fc Eggleton 



2001 



de Mink. Pols fc Hilditch 



20071 ). Therefore, the investigation of the evolution of massive bi- 



naries in case AR and AS is important to understand the formation of massive contact 
binaries, and then the formation of the rapidly rotating stars such as VFTS102. 

The purpose of this study is to investigate the binary merger channel for the progenitor 
of VFTS102 and to determine the detailed parameter range in which this channel produces 
VFTS102. Employing the Eggleton's stellar evolution code, we construct a grid of binary 
models for metallicity Z = 0.01 in Section 2 and 3, and then implement the results in a 
binary population synthesis study in Section 4. In Section 5, we give the discussion and 
conclusions. 



2 BINARY EVOLUTION CALCULATIONS 



In the binary merger channel for the progenitor of the fastest rotating 0-type star VFTS102, 
the primary of the initial detached binary would fill its Roche lobe and transfer some of its 
mass to the secondary. The secondary expands in response to the mass transfer from the 
primary. This system would evolve into contact if the secondary fills its own Roche lobe 
when two components are still MS stars. These contact systems would ultim ately merge 



into rotating stars due to the orbital angular momentum of the binary systems ( iJiang et al 



2010 



Langer 



2OI2I ). Furthermore, if both components are MS stars in case A binaries, the 



merge d stars would, be MS stars and evolve in a similar way to a normal star with that 
2OO9I : 



mass ( Chen fc Han 



Langerll2012l ). Therefore, the merged stars with the same mass as 



VFTS102 might be similar to VFTS102. 

There are some constraints on the progenitor of VFTS102 in the binary merger model. 
First, the detached binary, as the progenitor of VFTS102, should evolve into contact while 
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two component are still on the MS, which is largely determined by the initial period and mass 
ratio. Secondly, the total mass of the progenitor binary should be larger than or equal to 
the mass of VFTS102, which depends on the mass loss during t he evolution and th e merger 
process. The current mass of VFTS102 is approximately 25 (IDufton et al.ll201ll ). In this 
study, we do not consider the mass loss during the evolution and merger process and assume 
that the binary system as the progenitor of VFTS102 has total mass ~ 24.5 — 25.5 Mq. 

The merger process is complicated and the merger physics is still uncertain. Here, we 
adopt the following assumptions: (i) the binary systems merge immediately once both MS 
components fill their Roche lobes; (ii) the merged stars are homogeneously mixed; (iii) the 
system mass is conserved. The merged timescale (i.e. the time from a binary contact to 
merger) for massive contact binaries has remained unclear and we only adopt a simple 
assumption that the merger is instantaneous. As the components of contact binaries ro- 
tate rapidly during the merger process, the merged stars would be efficiently mixed. These 
merg ed stars migh t also undergo rotationally induced mixing during the s ubsequent evolu- 



2000 



tion jLangerll2012h. which is irnportant in the evolution o f massive stars (e.g. iHeger fc Langer 



iMaeder fc Meynet 



2000 



Howarth fc Smith 



20011 ). Therefore, it is reasonable to as- 



sume that the merged stars are homogeneously mixed. We roughly assume that the mass 
is conservative during the merger process. In fact, the binary systems might lose high an- 
gular momentum material duri ng the merge r process and the merged stars might rotate at 
subcritical rotation velocities (ILangerl |2012| ). The second and the third assumptions were 
often adopted in the study of the formation of blue stragglers by t he binary merger (e.g. 



Andronov. Pinsonneault fc Terndrup 



2006 



Chen fc Han 



2008 



20091) 



To determine whether the detached binary evolves into contact when two component 



are still MS stars, it is necessary to perform detailed binary evol ution calc u 



we use Eggleton's stellar evo 



ution code originally developed by lEggletonl (jl97ll . 119721 ): 



ations. Here 



Eggleton. Faulkner fcFlanneryl ( 119731 ). T his code has been updated with the 



physics during the last four decades fe.g. iHan. Podsiadlowski fc Eggleton 



1995 



1998 



Nelson fc Eggletonll200ll : lEggleton k. Kiseleva-Eggleton 



1994 



atest input 



Pols et al. 



2OO2I). We calculated our 



models with metallicity Z = 0.01 (hydrogen abundance X = 0.73 and Helium abundance 
Y = 0.26), which is close to the chemical composition of Large Magellanic Cloud {Z = 0.008). 



The mass transfer is determined to be dynamic when the mass-trans: 



lOM/txH? where ^kh is the thermal or Kelvin-Helmholtz timescale ( iNelson fc Eggleton 



er is greater than 
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200ll ). We construct a grid of stellar evolutionary models that covers the following ranges of 
initial primary mass Mio and initial mass ratio (go = ^^20/^^10): 

(1) 
(2) 
(3) 



logMio = 1.10, 1.15, ...,1.35, 
log(l/go) = 0.004, 0.05, 0.10, 0.95, 
log(Po/PzAMs) = 0.05, 0.10, ...,0.75, 



where Pzams is the period at whi ch the initially more mass ive component would just fill its 
Roche lobe on the zero-age MS (INelson fc EggletonI l200lh . The orbits are assumed to be 
circular. A binary is expected to become circularized during the RLOF by the tidal for ce on 
a timescale which is much smaller than the nuclear timescale ( Wang. Li &: Hani 2010 ). We 
assume that both total mass and orbital angular momentum are conserved in the evolution 
of the binaries. If the loss of the total mass and orbital angular momentum is considered, 
the orbit of the initial detached binaries would quickly shrink. More detached binaries with 
longer period would evolve into RLOF and contact, which could not evolve into contact in 
the conservation evolution. 



3 BINARY EVOLUTION RESULTS 

In Fig. 1, we present two representative examples of our binary evolution calculations that 
could evolve into contact while two components are still on the MS. It shows the evolutionary 
track in the radius-mass diagram, the mass-transfer rate and the evolution of the orbital 
period. Figs 1(a), (c) and (e) represent the evolution of a binary system with an initial mass of 
the primary of logMio = 1.20 (Afio =15.85 Mq), an initial mass ratio of log(l/go)=0.20 (go = 
0.63) and an initial orbital period of log(Po/-PzAMs) = 0.2 (Po = 1.498 d). The components 
start to evolve from point A and A', and the primary fills its Roche lobe on the MS at point 
B which results in case A RLOF (BC) as shown in Fig. 1(a). The secondary expands in 
response to the mass transfer (B'C) and fills its Roche lobe at point C. This system evolves 
into contact shortly after the start of RLOF. The mass-transfer rate from the primary 
to the secondary is not greater than lOM/txH as shown in Fig. 1(c). Therefore, the mass 
transfer is on a thermal t ime-scale and this case of binary evolution is classed as case AR 
(iNelson fc EggletonI l200lh . Fig. 1(e) shows that the orbital period of the system decreases 
from 1.498 to 1.37 d during the mass transfer phase because the mass is transferred from 
the massive component to the low-mass component. 



The progenitor of VFTS102 7 



0.9 



0.8 



0.7 



0.6 



M =15.85M 

10 o 

g„=0.63 

P =1.498day 



(a) 







1.55. 

1.50 

I 1-45. 
Q. 

1.40 
1.35 



12 
M{Mg) 



14 



16 



18 



(C) 



M =15.85M 

1 O 

C7„=0.63 

P =1 498day 



4 6 
age (Myr) 



(e) 



M =15.85M 

10 

<?o=0-63 

P =1.498 day 



4 6 
age (Myr) 



09 



08 

2. 07 
O 

0.6 



M =14.1 3M 

10 e 

q„=0.79 

P =1.142day 



1.5- 
1.4- 
?; 1.3- 

TO 
T3 

Q- 1.2- 

1.1- 




D\ (b) 



10 



12 
/W{Mg) 



14 



16 



18 



■ 1 1 1 1 1 1 1 1 1 1 1 

(d)- 






/H„=14.13Mg 






.7=0.79 

^0 






P„=1 .142day 







4 6 I 

age (Myr) 



10 



12 



■ 1 1 1 1 1 1 
M =14.13M 

10 o 


1 1 1 1 1 


90=079 




P^=1.142day 









10 



12 



age (Myr) 



Figure 1. Two examples of binary evolution calculations. In panels (a) and (b), solid curves show the evolutionary tracks 
in the radius-mass diagram for the primaries (ABC and ABCD) and the secondaries (A'B'C and A'B'C'D'), and the dotted 
curves represent the Roche lobe. The components start at point A or A', and evolve along curves until both components fills 
their Roche lobe at point C and C, or D and D'. In panels (c) and (d), the solid curves show the mass transfer rates and the 
dotted curves represent lOAf/tKH that is used to determined whether the mass transfer is dynamic. The evolution of orbital 
period is shown as solid curves in panels (e) and (f). 



Figs 1(b), (d) and (f) show another example for an initial mass of the primary of logMio 
= 1.15 (Mio =14.13M0), an initial mass ratio of log(l/go)=0.10 (go = 0.79) and an initial 
orbital period of log(Po/-PzAMs) = 0.1 (Pq = 1.142 d). The main difference between this ex- 
ample and the previous one is that the secondary fills its own Roche lobe, not during the rapid 
mass transfer phase (B'C), but during the subsequent slow mass trans fer phase (CD') as 
shown in Fig. 1(b). This case of binary evolution is identified as case AS (INelson fc Eggleton 



20011). In this C3iSG, clS shown in Fig. 1(d), the mass transfer rate is about 1.8 x 10 MQ/yr 



during the rapid mass transfer phase (BC) and about 5 x 10 ^Mo/yr during the slow mass 
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Figure 2. The outcomes of binary evolution calculations in the initial orbital period— mass ratio (log Po, log(l/qo)) plane while 
both components are still MS stars, where Po is the initial orbital period and qo is the initial mass ratio (for different initial 
primaries as indicated in each panel). Filled stars show tlic systems that evolve into contact during the rapid mass transfer 
phase while both components are on the MS (case AR). Open stars show the systems that evolve into contact during the slow 

mass transfer phase (case AS). Crosses indicate the systems that could not evolve into contact while two components are still 
on MS. Solid lines show the boundaries (without left boundaries) for the initial parameters for which might be the progenitor 
of VFTS102. 



transfer phase (CD). The orbital period of this system first decreases from 1.142 to 1.097d, 
and then increases to 1.478 d after the mass ratio is reversed as shown in Fig. 1(f). 

Fig. 2 summarizes the outcomes of binary evolution calculations in the initial orbital 
period— mass ratio (logPo? log(l/5'o)) pl^^iic while both components are still MS stars. Filled 
stars show the systems that evolve into contact during the rapid mass transfer phase while 
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both components are on the MS (case AR). Open stars show the systems that evolve into 
contact during the slow mass transfer phase (case AS). Crosses indicate the systems that 
could not evolve into contact while two components are still on MS. We also present the 
boundaries (solid lines) for the initial parameters for which might be the progenitor of 
VFTS102. The right boundaries are set by the condition that the system could evolve into 
contact while two components are still on the MS. The upper and lower boundaries are 
caused by the constraint of the mass of VFTS102. The left boundaries are not plotted in 
Fig. 2, which are set by the period (Pzams) at which the initially more massive component 
would just fill its Roche lobe on the zero-age MS. It is seen from Fig. 2 that the progenitor 
binary of VFTS102 with go ^0.7 (log go ^0.2) should have an initial orbital period shorter 
than 3.76 — 4.25 days and they would evolve into contact during the rapid mass transfer 
phase. The progenitor binary of VFTS102 with go > 0.7 (log go <0.2) should have an initial 
orbital period shorter than 1.44 — 1.55 days and they would evolve into contact during the 
rapid or slow mass transfer phase. 



4 BINARY POPULATION SYNTHESIS 



We have performed a series of detailed Monte Carlo simulations to investigate the binary 
merger channel. In each simulation, we follow the evolution of 100 million sample binaries 
(very wide binaries are actually single stars) according to grids of stellar models of metallicity 
Z = 0.01 and the evolution channels described above. We adopt the following input for 
the simulations: the initial mass function (IMF) of the primaries, the initial mas s -ratio 



distribution, and the distri 



Lii] 



2009 



Wang. Li fc Han 



Dution of initial orbital separations (e.g. 



Han et al. 



12002 



2003 



2010): 



(i) We use a simple approximation to the IMF of 



Miller &: Scalo (1979) and the mass of 



the primary is generated using a formula of lEggleton. Tout &: FitechettI (1l989[ l 
0.19X 



M- 



10 



(1 - X)0-75 + 0.032(1 - X)0-25 ' 
where X is a random number u niformly distributed between and 1. The study of IMF by 
Kroupa. Tout fc Gilmord ( 119931 ) supports this IMF. 
(ii) We take a uniform mass-ratio distribution (Set 1), 



n(go) = 1, < go < 1 



(5) 
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Table 1. The results of the binary population synthesis for the binary merger channel. 



Set 


i(<?o) 


iVm 




/b 






1 


Uniform 


1006 


11581 


8.7% 


4875 


17.1% 


2 


Rising 


787 


10834 


7.3% 


4223 


15.7% 


3 


Uncorrelated 


52 


7559 


0.7% 


4830 


1.1% 



Columns: Set-name of different simulation set; n(go)-the distribution of the mass ratio distribution; A^m-the number of binaries 
with total mass ~ 24.5 — 25.5 Mq that would merge into single stars; JV^-the number of the initial binaries with total mass 
~ 24.5 — 25.5 Mq; /ij-the fraction of binaries that would merge into single stars (f^ = N^/N^)\ Afg-the number of the initial 
single stars with mass ~ 24.5 — 25.5 Mq; /a-the fraction of single stars that might be produced from the binary merger channel 

(/s = N^/(N^ + N,)). 



where go = M20/M10 flMazeh et alJll992l : IColdberg fc Mazehl Il994h . In order to study the 



influence of the mass ratio distribution, we also take a rising mass ratio distribution (Set 2) 
n{qo) = 2qo, < go < 1 (6) 

and an alternative mass-ratio distribution where both components are chosen randomly and 
independently from the same IMF (Set 3). 

(iii) We assume that all stars are members of binary systems and that the distribution of 
separations is constant in log a for wide binaries, where a is the orbital separation and falls 
of smoothly at small separation 

-plao^ ' - °' (7) 

Osep, ao < a < ai, 

where ctsep ~ 0.070, oq = lOR©, ai = 5.75 x 10^ R© = 0.13 pc, and m ^ 1.2. This distribution 
implies that the numbers of wide binary systems per logarithmic interval are equal, and that 



abou t 50% of stellar systems have orbital periods less than 100 yr (IHan. Podsiadlowski fc Eggleton 



19951 ). A circular orbit is assumed for all binaries. 



The results of the binary population synthesis are shown in Table 1 for three sets of 
simulations. A/'m, N\, and As are the number of binaries that would merge into single stars, 
the initial binaries and the initial single stars, all of which are in the range of mass (or total 
mass for the binaries) from 24.5 to 25.5 M©. /b is the fraction of binaries that would merge 
into single stars, and fs is the fraction of single stars that might be produced from the binary 
merger channel, where /b = N^/Ny,, and fs = A'm/(Am + A's)- Set 1 is our standard model 
for the binary merger channel and we vary the mass ratio distribution in the other sets to 
examine its influence on the final results. The simulation for the binary merger channel gives 
/b ~ 8.7% and fs ~ 17.1% according to our standard model (Set 1). This suggests that about 
8.7% of binaries would evolve into contact and merge into rapidly rotating single stars, and 
about 17.1% of single stars might be produced from the binary merger channel and should 
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0.0 0.2 0.4 0.6 0.8 1.0 



log (days) 

Figure 3. The relative number distribution of the progenitor of VFTS102 in the initial orbital period— mass ratio (logPo, 
log(l/go)) plane. 

have similar properties to VFTS102. The simulation for a rising mass-ratio distribution (Set 
2) gives /b ~ 7.3% and /g ~ 15.7%, which are close to the results in Set 1. If we adopt a 
mass-ratio distribution with uncorrelated binary components, two fractions (/b ~ 0.7% and 
/s ~ 1.1%) are much lower than the results in Set 1. This is because most of binary systems 
have too small a mass ratio (go) to locate in the region of the progenitor of VFTS102 for 
different initial primaries. 

Fig. 3 shows the relative number distribution of the progenitor of VFTS102 in the initial 
orbital period— mass ratio plane (Set 1). It is found that the distribution of the progenitor 
of VFTS102 increases with decreasing orbital period and increasing mass ratio (go)- This 
suggests that VFTS102 mainly comes from the binary with short period and large mass 
ratio. The most likely progenitor binary of VFTS102 might have an orbital period Pq < 1.5 d 
and a mass ratio go > 0.63, which roughly corresponds to a binary with a primary mass 
(Mio ~ 12.5 — 15 M0) according to Fig. 2. This is mainly because the number of initial 
binaries increases with the decreasing primary mass according to the IMF of the primaries, 
and more systems with low-mass primaries could be located in the production region of 
VFTS102. 

5 DISCUSSION AND CONCLUSIONS 

In this paper, we investigated the binary merger channel for the progenitor of the fastest 
rotating 0-type star VFTS102. We carried out detailed binary evolution calculations and 
obtained the initial parameters of the progenitor of VFTS102 in the orbital period- mass 
ratio plane. By performing Monte Carlo simulations, we estimated the fraction of binaries 



12 Jiang et al. 

that would merge into single stars and the fraction of single stars that might be produced 
fro m the binary merger channel. 



Dufton et al 



(1201 ll ) proposed that VFTS102 originated from a binary system, and the 
mass transfer from its component resulted in the fast rotation of VFTS102. They sug- 
gested that the subsequent supernova explosion of its component kicked both components 
and led to an anomalous radial velocity for VFTS102. We showed that the fast rotation 
of VFTS102 could be understood within the framework of the binary merger channel. In 
this channel, a detached binary evolves into contact and then merges into a r apidly ro- 



tating single star due t o the orbital angular mornentum of the binar y systems (jWebbink 



1976 



Jiang et al. 



2010 



de Mink. Langer fc Izzard 



2011 



Langeiil2012l ). The anomalous ra- 



dial velocity for VFTS102 might be produce d by the progenitor binary, which might be 



a binary dynamically ejected from a cluster (Ide Wit et al. 



2005 



2011 



Gvaramadze fc Gualandris 



2011 



Eldridge. Langer fc Tout 



)• 



We mapped out the initial parameters of the proge nitor of VF TS102 in the binary merger 
channel. Our results are very similar to those given by iPold (1l994l ). He showed the parameter 
space of initial mass ratio and orbital period for the systems with 16 Mq, which was explored 
for the formation of contact binaries. By performing Monte Carlo simulations, we found that 
the most likely progenitor of VFTS102 might have a primary mass Mio ~ 12.5 — 15 Mq, a 
mass ratio go > 0.63 and an orbital period Pn < 1.5 d. 



According to the work of iPodsiadlowski. Joss fc Hsd (119921 ). the fraction of binaries that 



experience case A and merge is about 14.8%, which is larger than 8.7% found in this study. 
This is because they assumed that case A mass transfer always leads to a merger of two 
component, in which some of binaries would reach c ontact while one or both components 



evolve past the terminal MSjjNe 



son &: Eggletonll200ll ). and these binaries are not considered 



in our study. In addition. iLangerl (120121 ) suggested that up to 10% or more of the stars would 
be merger remnants, which is close to the result of our study (17.1%). Therefore, the binary 
merger channel might be one of the main channels for the production of rapidly rotating 
star like VFTS102 and it is necessary to consider this channel for the study of the formation 
of rapidly rotating stars . 



Tylenda et al 



(I2OIII ) showed the direct observation evidence that contact binaries in- 
deed merge into a single object from the identification of V1309 Sco. The observations 
also show some possible progenitor of VFTS102-like objects that are short-period , massive 
semi-detached or contact binaries, such as V Pup (P=1.45 d lAndersen et al.lll983l ). lU Aur 
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fP=1.811d 



Harries. Hilditch fc Hil]||l998l). V606 Cen ('P=1.495 diLorenz. Maver. fc Drechsel 



19991) and TU Mus (P=1.387d 



Terrell et al. 



2003 



Qian et al.ll2007l ). In addition, some mass 



should be lost fror n binary systems during the na erger process to carry away the excess an- 



gular momentum (jjiang et al. 



2010 



Langerll2012l ). which might form circumstellar material 



or a circumstellar disc. This cou l d exp lain the observation that VFTS102 has circumstellar 



material found by 



Dufton et al. 



(1201 If ). 

In this paper, we investigated the binary merger channel for the progenitor of VFTS102 
without considering the mass loss during the evolution and merger phase. However, the mass 
loss is an important parameter in the evolution, the mass transfer phase and the merger phase 
of massive binaries, which needs to be considered in the future study. In addition, the merger 
process of the binaries and the evolution of the merger stars are still uncertain and need to 
be investigated further. 
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